Dissipative dynamics of superconducting hybrid 
qubit systems 

Enrique Montes, Jesiis M Calero and John H Reina 

Departamento dc Fi'sica, Universidad del Valle, A. A. 25360, Cali, Colombia 
E-mail: enriquem@univalle.edu.co and jhreina@univalle.edu.co 

Abstract. We perform a theoretical study of composite superconducting qubit systems for 
the case of a coupled qubit configuration based on a hybrid qubit circuit made of both charge 
and phase qubits, which are coupled via a, <^ ctz interaction. We compute the system's 
eigen-energies in terms of the qubit transition frequencies and the strength of the inter-qubit 
coupling, and describe the sensitivity of the energy crossing/anti-crossing features to such 
coupling. We compute the hybrid system's dissipative dynamics for the cases of i) collective and 
ii) independent decoherence, whereby the system interacts with one common and two different 
baths of harmonic oscillators, respectively. The calculations have been performed within the 
Bloch-Redfield formalism and we report the solutions for the populations and the coherences 
of the system's reduced density matrix. The dephasing and relaxation rates are explicitly 
calculated as a function of the heat bath temperature. 
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1. Introduction 

A quantum computer uses quantum two-level systems or quantum bits (qubits) which must 
be essentially interacting systems weakly or, even more desirable, completely decoupled from 
dissipative environments [Ij. The low temperature superconductors, like aluminum, are good 
candidates for this purpose. For this reason, in the last few years, interest in the study of 
superconducting-qubit systems has grown O [3]. In these kind of systems, the circuit must be 
cooled to a temperature at which the thermal energy = ksT is less than the transition 
energy hujo^i between the qubit basis states, say |0) and Another important requirement 
is that hujQ^i << A, where A is the energy gap of the superconducting material [3]. The basic 
elements used in the preparation of superconducting qubits are the Josephson junctions, which 
are one of the few electronic systems which, besides being nonlinear, also exhibit the property of a 
very weak dissipation [5] . Quantum information processors require a qubit coherent dynamics in 
order to allow the construction of quantum logic gates and circuits [Ij. The hybrid systems have 
had an impact in the development of these circuits because they can bind the scales of atomic 
motion to macroscopic degrees of freedom [6]. Some solid-state systems, such as the so-called 
superconducting qubits, allow the implementation of such hybrid systems in their respective 
circuits which can be connected in different ways, in the form of simple circuit elements [7|- 

2. Model 

Hybrid quantum circuit. — A generic superconducting qubit can be modeled by the Hamiltonian 

Hqb = -^^^z - ^'^^x, (1) 




where e and A denote the longitudinal and 
transversal parameters of the correspond- 
ing qubit, and fTj are the usual Pauli ma- 
trices. Charge and flux qubits are two dif- 
ferent types of superconducting qubits for 
quantum computing. In short, the charge 
qubit has the advantage of a more flexible 
controllability via external parameters: it 
can be conveniently controlled by a voltage 
gate or an applied magnetic flux. These ex- 
ternal control parameters appear in the lon- 
gitudinal (az) and transverse (ax) terms of 
the circuit's reduced Hamiltonian. For the 
flux qubit, the longitudinal term can be con- 
trolled by the applied magnetic flux, but it 
is hard to control the transversal term via 
an external parameter. 

The hybrid quantum circuit we study in 
this work is a system that couples a charge and a flux qubit, as schematically shown in Figure 
[H with an effective interaction due to a Josephson junction that binds them. This system has 
been proposed in order to control the transversal term of the flux qubit with the charge qubit 
[7]. A calculation of the explicit Hamiltonian for the hybrid qubit system reads 



Figure 1. Hybrid quantum circuit: charge (left) 
and flux (right) qubits are effectively coupled due 
to the Josephson junction Ejq. Eji denotes the 
Josephson energy of each junction (crossed boxes). 
The charge (i = 1) and the flux (i = 2) qubits are 
crossed by externally controlled magnetic fluxs 
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where ei and Ai are the longitudinal and transversal parameters of the charge qubit, and likewise 
62 and A2 for the flux qubit. The interaction between the qubits gives rise to an effective ax^cTz 
geometric term [7] with strength x- 



3. Coupling to the bath 

In order to study the system's dissipative dynamics, we modeled the dissipative environment 
as either one or two baths of harmonic oscillators which are coupled to the cj^ components of 
each qubit. The first case considers each qubit coupled to its own bath of oscillators in the form 
a^z^X^, where X* [i = 1, 2) denotes the bath coordinates. The full interacting Hamiltonian reads 

H2B = E {-^^^^^ - l^^^^^ + l^^^^') + X'y^^cjf) + Hbi + Hb2, (3) 

where Hsi {i = 1,2) is the Hamiltonian of the i-th bath. In the second case, we consider that 
the qubits are coupled to a common bath. Hence, the full system's Hamiltonian reads 




where X is the bath coordinate and Hb represents the Hamiltonian of the common bath. 

Block- Redfield formalism. — This formalism is derived from a projector operator approach and 
provides an important tool for finding a set of coupled master equations which describe the 
dynamics of the reduced density matrix for a given system in contact with a dissipative environ- 
ment. The Liouville equation for the total density operator px (of the whole system) is given by 



ih^PT{t) = [H{t), pT{t)]. The corresponding reduced density matrix obeys, in the Born-Markov 
hmit, the so-called Redfield equation [H]: 

P = —i^nmPnm{t) — ^ RnmklPkl{t), (5) 
kl 

where the first term represents the reversible motion in terms of the transition energies 
hiOnm = En— Em, and the second term describes the system's relaxation. The Redfield relaxation 
tensor Rnmki includes the dissipative efects of the system-environment coupling: 

Rnmkl = ^Im X] ^nrlk + ^nk X] '^Lrm. ~ '^\mnk ~ ^\mnk^ (^) 
r r 

and its components are given by the 'golden rule' expressions 

^\tlk = r dte-'^-^'{Hi,Ut)Hi^nk{^)), 

J 

rJ;;L = n-'J^^dte-'^^'-'{Hj,UO)Hi,nkit))- (7) 

Here, Hj{t) = exp{ ^^jP'' )Hj exp( ~^^^* ) is the interaction term in the interaction picture, and 
the bracket denotes the thermal average over the bath's degrees of freedom. 



4. Results and Discussion 

System's eig en- energies. — These are worked out in the singlet /triplet basis states. In this basis, 
the Hamiltonian H2q of the two-qubit system assumes the matrix form: 
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where we have introduced e = ei + £2, Ae 



qubits with the same transition energies ei = 62 and transversal coupling Ai 
Ae = Ar/ = 0. 

In Figure [21 we present the system's energies for three different scenarios, in which two of the 
parameters e, rj, x have been fixed and the third one is varied. In Figure [2)[|a) we have varied the 
longitudinal parameter e, showing a crossing between the energy of the singlet {Ei) and one of 
the energies of the triplet (E^). In Figure [2]|^b) the infiuence of the transversal parameter rj has 
been plotted. We notice the appearance of the same crossing observed in Figure [2l||a) , between 
the singlet and the £^3-triplet state, which indicates that these states do not interact. Figure 
W^c) shows the dependence on the coupling term x- Now, a crossing between the energy of the 
singlet state and that of the triplet state E4 takes place. Note that for large values of x, there 
is a degeneracy between the singlet state Ei and the triplet state E^. This feature means that 
for large positive values of the parameter x, a large ferromagnetic coupling is present between 
the qubits. 
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Dissipative dynamics and temperature effects. — The solution of the Redfield equation pnm = 
— J2ki RnmkiPkiit) allows to analyse the system's reduced dynamics, for which explicit 
expressions have been found for the components of Equation ([7]), for both types of coupling 
to the environment. In the case of independent coupling we get 
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Figure 2. Energy spectrum of the coupled hybrid two-qubit system, (a) x = 10 GHz, rj = 1 
GHz, and e is varied, (b) x = 1 GHz, e = 1 GHz (r/ is varied), and (c) e = r/ = 1 GHz (x is 
varied). In all the graphs, Ei is the singlet state and E2, E^, and E^ are the triplet states. 



where the a^'s correspond to the bath coupling parameters for the specific cases of Ohmic 
spectral densities Ji{u>) = aihio/(^l + and a^^ij = {Ei\az \Ej). For the case of collective 

coupling, J{io) = ahLo/(^l + and the rates become 
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Figure 3. Populations dynamics, Pii{t). The density matrix calculation has been performed for 
the initial conditions: a) pure state and b) maximally entangled Bell state l^*"*"). Lower 

graphs show the relaxation rates (see main text) as a function of the bath temperature for the 
cases of c) independent and d) collective coupling to the environment. Ui = a = 10~^. 



The reduced density matrix has been computed from the system of differential equations 
([5]). Figure [3] shows the time dependence of the system's populations pa for a coupling to two 




different kinds of environments: i) independent and ii) collective or common baths. Figures 3(a) 
and 3(b) correspond to the initial conditions \Ei), and |<I>+) = (|00) + |11)), respectively. 
The population pn has a slightly faster decay in the case of collective coupling for the initial 
condition \Ei), whereas this behaviour is reversed for the case of the entangled initial condition 
Thus, the dissipative population dynamics is very sensitive to initial state preparations. 
The behaviour of the relaxation rates Vll^ = Re{Rnnnn) as a function of the bath temperature 
is shown in Figures [3||c) and (d), for the initial conditions \Ei). The superscripts i, and c stand 
for independent (Fig. E^c)) and collective (Fig. E^d)) couplings respectively. In agreement with 
graph 3(a), this shows that the rate Fn tends to a higher value for the case of collective coupling. 

In FigureHl we have plotted the dephasing (decoherence) rates Fn,m = Re{Rnmnm), as a function 
of the bath temperature, for the cases of a) independent and b) collective coupling. We note 
that, in both cases, the decoherence rates are of the same order of magnitude. Furthermore, it 
is observed that is the minimum dephasing rate, while F31 displays the higher rate. 

By comparing the decoherence and the relaxation rates obtained, we conclude that both 
types of dissipation mechanisms occur on the same time scale. This implies a long decoherence 
time, which is a key requirement in the construction of quantum devices that are expected to 
work as quantum information processors. 
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